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During  the  period  of  the  grant  we  continued  our  modeling  and  analysis  of 
the  mechanical  behavior  of  complex  composite  materials.  We  have  developed 
analytical  and  numerical  modeling  techniques  of  the  influence  of  general  com¬ 
posite  laminate  orientation  on  the  ultrasonic  behavior  of  anisotropic  plates 
and  substrates.  In  the  second  phase  we  introduced  piezoelectric  effects  into 
our  modeling.  In  the  final  phase  we  studied  the  dynamic  response  of  the  lay¬ 
ered  composite  to  transient  loadings  in  the  form  of  time  dependent,  source 
loads.  Complete  documentation  of  the  results  of  the  first  two  phases  were 
submitted  to  the  AFOSR  in  the  form  of  yearly  reports.  Since  these  reports 
contained  completed  items  which  also  appeared  in  the  literature  as  archival 
publications,  we  need  not  rereport  them  here.  Results  for  the  recently  com¬ 
pleted  works  have  also  been  prepared  for  journal  publication.  Description  of 
these  works  are  included  here  in  details. 


Chapter  1 


WAVE  PROPAGATION  IN 
ANISOTROPIC  MEDIA  DUE 
TO  INTERNAL  HARMONIC 
LINE  LOADS 


1.1  Introduction 

Plane  harmonic  wave  interaction  with  homogeneous  elastic  anisotropic  media, 
in  general,  and  with  layered  anisotropic  media,  in  particular,  have  been  ex¬ 
tensively  investigated  in  the  past  decade  or  so.  This  advancement  has  been 
prompted,  at  least  from  a  mechanics  point  of  view,  by  the  increased  use  of 
advanced  composite  materials  in  many  structural  applications.  Being  both 
anisotropic  and  dispersive,  composite  materials  required  indepth  understand¬ 
ing  of  their  mechanical  behavior.  The  list  of  relevant  literature  is  rather  long 
and  we  thus  refer  the  reader  to  selective  recent  works  for  further  representative 
references  [1-6].  Some  of  the  difficulties  inherent  in  the  treatment  of  wave  prop¬ 
agation  in  anisotropic  media  can  be  illustrated  by  considering  slowness  wave 
surface  methods  for  infinite  media  [1,7-9]  (Synge.  1957;  Fedorov.  1968:  Mus- 
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grave,  1970).  According  to  this  analysis,  the  slowness  surface  of  an  isotropic 
material  consists  of  two  concentric  spherical  sheets,  an  inner  one  representing 
the  longitudinal  wave  mode  and  an  outer  one  representing  the  two  degenerate 
transverse  modes.  In  the  anisotropic  case,  we  find  there  are  three  general  wave 
surfaces,  one  for  the  so-called  quasi-longitudinal  wave  and  two  nondegenerate 
sheets  for  the  quasi-transverse  waves.  Moreover,  the  surfaces  will  no  longer 
be  spherical  in  shape  but  will  reflect  the  elastic  symmetry  or  assy  met  ry  of  the 
material.  Perhaps  the  most  severe  consequence  of  elastic  anisotropy  in  infinite 
media  is  the  loss  of  pure  wave  modes  for  general  propagation  directions.  This 
fact  also  implies  that  the  direction  of  energy  flow  |i.e.,  group  velocity)  does 
not  generally  coincide  with  the  wave  vector,  or  wavefront  normal.  Clearly, 
uncoupled  pure  potentials  (such  as  are  found  in  the  isotropic  case)  are  much 
simpler  to  treat  than  the  mixed  modes  characteristic  of  anisotropic  materials. 
For  wave  propagation  in  directions  of  symmetry  some  wave  types  revert  to 
pure  modes,  leading  to  a  simpler  characteristic  equation  of  lower  order.  A 
key  condition  which  we  found  to  facilitate  our  previous  analysis  is  that  at 
any  boundary  all  wave  vectors  must  lie  in  the  same  plane.  This  requirement 
implies  that  the  response  of  the  media  will  be  independent  of  the  in-plane 
coordinate  normal  to  the  propagation  direction.  We  therefore  conducted  all  of 
our  analysis  in  a  coordinate  system  formed  by  the  line  load  direction  and  its 
normals  rather  than  one  determined  by  material  symmetry  axes.  This  choice 
lead  to  a  significant  simplification  in  our  algebraic  analysis  and  computations. 
Compared  to  the  extensive  literature  on  the  interaction  of  plane  harmonic 
waves  with  anisotropic  media,  very  little  work  is  available  on  the  response  of 
such  media  to  concentrated  source  loading.  Here  concentrated  sources  include 
point  as  well  as  line  loads.  Harmonically  pulsating  and  transient  sources  are 
common  types  of  such  loading.  Understanding  the  response  of  elastic  solids  to 
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internal  mechanical  sources  lias  long  been  of  interest  to  researchers  in  classi¬ 
cal  fields  such  as  acoustics,  seismology,  as  well  as  modern  fields  of  application 
like  ultrasonics  and  acoustic  emission.  It  is  known  that  whenever  a  material 
undergoes  a  local  failure,  elastic  waves  are  generated  due  to  the  rapid  release 
of  localized  strain  energy.  Such  radiation,  for  example,  is  known  as  acoustic 
emission  in  the  field  of  nondestructive  testing  of  materials.  In  seismology  it  is 
of  course  known  by  the  earthquake. 

A  quick  review  of  available  literature  on  this  subject  reveals  that  most  of 
the  work  done  so  far  is  carried  out  on  isotropic  media.  The  effect  of  imposed 
line  load  in  homogeneous  isotropic  media  has  been  discussed  by  several  in¬ 
vestigators  ever  since  Lord  Rayleigh  discovered  the  existence  of  surface  waves 
on  the  surfaces  of  solids  [10].  An  account  of  the  literature  dealing  with  this 
problem  through  1957  can  be  found  in  Ewing,  Jardetsky  and  Press  [18].  Most 
of  the  earlier  work  [13-15]  followed  Lamb  [11,12].  who  apparently  was  the  first 
to  consider  the  motion  of  semi-infinite  space  caused  by  a  vertically  applied 
line  load  on  the  free  surface  or  within  the  medium.  He  was  able  to  show  that 
displacements  at  large  distance  consists  of  a  series  of  events  which  corresponds 
to  the  arrival  of  P-,  S-,  and  Rayleigh  surface  wave.  The  analytical  approach 
used  in  the  above  mentioned  investigations  and  others  [13-15]  can  be  summa¬ 
rized  as  follows  ;  the  steady  state  problem  for  harmonic  waves  propagating 
in  infinite  isotropic  media  is  solved  at  first  and  then  generalized  to  the  case 
of  half-space  using  superposition  technique.  For  transient  source  loading  re¬ 
sults  can  be  obtained  from  those  corresponding  to  harmonic  ones  by  a  Fourier 
integral  approach. 

In  this  paper,  we  closely  follow  the  formal  developments  in  our  previous 
works  [1-6]  and  study  the  response  of  several  anisotropic  systems  to  harmoni¬ 
cally  pulsating  buried  line  loads.  These  include  infinite,  semi-spaces  and  plate 
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systems.  Our  analysis  will  be  carried  out  on  anisotropic  media  possessing 
monoclinic  or  higher  symmetry.  The  load  will  be  in  the  form  of  a  normal 
stress  load  acting  at  an  arbitrary  direction  within  the  materials  in  the  plane 
of  symmetry.  We  use  a  building  block  approach  in  which  we  start  by  driving 
results  for  an  infinite  media.  Subsequently  we  obtain  the  results  for  half-spaces 
by  using  superposition  of  the  infinite  medium  solution  together  a  scattered  so¬ 
lution  from  the  boundary.  The  sum  of  both  solutions  has  to  satisfy  stress  free 
boundary  conditions  thereby  leading  to  complete  solution.  Lastly  we  proceed 
to  develop  solutions  for  a  plate  by  insuring  both  the  infinite  media  solution  and 
the  scattering  solution  to  satisfy  the  two  free  surface  of  the  plate.  Our  work 
will  be  easily  executed  by  using  the  linear  transformation  approach  in  which 
we  identify  the  line  load  with  the  x2  direction.  This  implies  that  all  involved 
field  variables  will  be  independent  of  the  x2  direction.  Nevertheless,  and  in 
general  we  will  have  three  nonvanishing  particle  displacements.  Material  sys¬ 
tems  of  higher  symmetry,  such  as  orthotropic,  transversely  isotropic,  cubic, 
and  isotropic  are  contained  implicitly  in  our  analysis.  The  equivalent  crystal 
systems  of  monoclinic,  orthohombic,  hexagonal  and  cubic  may  be  substituted 
for  the  elastic  material  systems  analyzed  here.  We  demonstrate  numerical  re¬ 
sults  results  drawn  from  concrete  examples  of  materials  belonging  to  several 
of  these  symmetry  groups.  For  orthotropic  and  higher  symmetry  materials 
where  the  remaining  two  principal  axes  lie  in  the  plane  of  the  plate,  the  parti¬ 
cle  motions  in  the  saggital  and  the  normal  to  it  uncouple  if  propagation  occurs 
along  either  of  these  in-plane  axes. 

1.2  Theoretical  development 

Consider  an  infinite  anisotropic  elastic  medium  possessing  monoclinic  symme¬ 
try.  The  medium  is  oriented  with  respect  to  the  reference  cartesian  coordinate 
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system  =  (x^.Xj.xij)  sudi  that  the  x'3  is  assumed  normal  to  its  plane  of 
symmetry  as  shown  in  figure  1.  The  plane  of  symmetry  defining  the  mono¬ 
clinic  symmetry  is  thus  coincident  with  the  x\  —x'2  plane.  With  respect  to  this 
primed  coordinate  system,  the  equations  of  motion  in  the  medium  are  given 

by  [i] 


and,  from  the  general  constitutive  relations  for  anisotropic  media. 


® ij  Cijkle’kli 


i,j,k,l  =  1,2,3 


by  the  specialized  expanded  matrix  form  to  monoclinic  media 


where  we  used  the  standard  contracted  subscript  notations  1  — »  11,2 


(1.3) 

_>  22, 


3  — »  33,  4  — *  23,  5  — >  13,  and  6  — ►  12,  to  replace  Cijhi(i,j,k,l  =  1.2,3)  with 
cPq(p,q  =  1,2,. ..,6).  Thus,  c45  stands  for  c23i3,  for  example.  Here  cr'x] , e';  and 
u\  are  the  components  of  stress,  strain  and  displacement,  respectively,  and 
p'  is  the  material  density.  In  Eq.  (3),  7 =  2etJ  (with  i  ^  j  )  defines  the 
engineering  shear  strain  components. 

In  what  follows,  we  study  the  response  of  the  infinite  medium  to  a  uniform 
harmonic  line  load  applied  along  a  direction  that  makes  an  arbitrary  azimuthal 
angle  <f>  with  the  x\  axis.  That  is,  the  direction  <j>  =  0  coincides  with  the 
reference  coordinate  x\.  Since,  as  was  pointed  out  in  the  introduction,  the 
response  of  the  medium  to  such  a  wave  is  independent  of  the  applied  line 
direction,  we  conduct  our  analysis  in  a  transformed  coordinate  system  x,  = 
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(xj  ,£2,  £3)  formed  by  a  rotation  of  the  plane  x't  —  x2  through  the  angle  o 
about  the  X3  direction.  Thus,  the  direction  x2  will  coincide  with  the  line  load 
direction. 

Since  c,jkt  is  a  fourth  order  tensor,  then  for  any  orthogonal  transformation 
of  the  primed  to  the  non-primed  coordinates,  i.e.,  it  transforms  according  to 

C\jkl  ^im  $ jnfiko&lpt'rrtnop  (1-d) 


where  is  the  cosine  of  the  angle  between  x'  and  x,  ,  respectively.  For  a 
rotation  of  angle  o  in  the  ij  —  x2  plane,  the  transformation  tensor  reduces 
to 

cos  0  sin  <j>  0 

0ii  —  -sin  cj)  cos  <p  0  (1.5) 

0  0  1 

If  the  transformation  (5)  is  applied  to  Eq.(3),  one  gets 


Cu 

C12 

Cl3 

0 

0 

C\2 

C-2  2 

C23 

0 

0 

Cl3 

C23 

C33 

0 

0 

0 

0 

0 

C44 

c45 

0 

0 

0 

C45 

c55 

C16 

c26 

C36 

0 

0 

(1-6) 


where  the  relations  between  the  cvq  and  c'pq  entries  are  listed  in  [5j.  Notice 
that,  no  matter  what  rotational  angle  <f>  is  used,  the  zero  entries  in  Eq.(3)  will 
remain  zero  in  Eq.(6).  In  terms  of  the  rotated  coordinate  system  x/j,  we  also 
write  the  momentum  equations  as 


d<Tjj  .  _  82u , 
dxj  1  ^  dt 2 


As  mentioned  earlier,  in  the  rotated  system,  the  elastic  wave  equations,  for 


wave  propagating  in  the  Xi  —  x3  plane,  are  independent  of  x2  .  Nevertheless, 


the  particle  motion  can  generally  have  three  nonzero  components  and 
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u3.  We  identify  the  u>  displacement  as  belonging  to  th*.  horizontally  polarized 
(SH)  wave.  Here,  the  equations  describing  these  three  wave  motions  are  cou¬ 
pled.  complicating  the  analysis  of  free  waves  in  the  medium.  In  terms  of  the 
displacement  components  the  equations  are  written  in  the  expanded  form 


,  d2  d2  d2  ()  2  . 

\-Cn~d72  +  C55CTUl  +  ^Cui7T, 1  +  C457T^jU2 


'OrV 


'  dx2 


d  f  0 

+  ^ - [( c13  +  C55 )  — - j  U3 

OX3  US  1 

02  d2  d2  ()2 

Ic,657f  +  +  +  C44&y“; 

a  „  i)  , 

+d7r,lici‘  +  c,')o7l ‘1U3 

d  r,  .  d  d  r/  c) 

^((C13  +  Css)^lu‘ +  a^l(C36  +  c<s  W“2 

r  ^  1 

+|CM&| + 


d2Uy 

PHF  ~  J'  (Sa) 


d2V"2 

p-^r-hm 


with 

fi  -  Qi6{xl)6(x3)F(t)  ( 1.9) 

1.3  Source  characterization 


Using  the  above  geometric  arrangement  we  proceed  to  describe  the  line  load 
by  a  direc  delta  body  force  function  located  at  the  arbitrary  location  (x3  =  x3 
).  This  choice  of  source  location  is  arbitrary  and  is  not  necessary  for  the  cases 
of  propagation  in  either  the  infinite  nor  the  semi-infinite  media.  However,  as 
will  be  shown,  this  choice  will  lead  to  great  simplification  in  the  algebraic  ma¬ 
nipulation  for  the  propagation  in  the  plate.  Solutions  to  the  present  problem 
can  be  accomplished  by  following  the  procedure  used  by  Achenbach  [20].  Ac¬ 
cording  to  this  procedure,  the  infinite  space  can  be  thought,  of  as  consisting 
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of  two  semi-spaces  whose  artificial  interface  contains  the  applied  load.  Thus 
the  upper  half-space  occupies  the  region  x3  >  .r3  whereas  the  lower  half-space 
occupies  the  region  x3  <  As  a  result  of  this,  appropriate  condition  must 
be  specified  at  the  artificial  interface. 

In  order  to  be  able  to  specify  these  int  rface  conditions  we  consider  a  very 
thin  "interface"  layer  extending  from  x3  =  jHIJ  —  o  to  ,r3  =  4-  o.  We  then 

start  by  requiring  continuity  on  the  displacement  components  and  their  time 
derivatives,  namely 


u,(xt.  j3.0|^ 
fj(n)u,(ji..r3.Q  zp+0 


=  o 
=  o. 


i  =  1.2,3 


( 10a) 


(10b) 


Integrating  Eqs.(8a,b.c)  across  the  interface  layer  and  using  the  continuity 
equations!  10a. b),  leads  to  the  discontinuity  in  the  displacement  spatial  deriva¬ 
tives  resulting  in  the  following  jump  conditions. 


du, 

(CSSTT— 

C/CC3 

,  dvl\\*P3+° 

+c«te)l4-.  = 

-QiS{xx)F{t) 

(1U) 

du  1 

(C4S&, 

+*«jr>i£:  = 

0x3  *  0 

-Q26MFV) 

(lib) 

du3  xP  u0 

-Q3S(x,)F(t) 

(He) 

Now,  let’s  consider  the  special  case  where  Q\  =  Q 2  =  0  and  Q3  —  Q,  For  this 
situation  equations  (11)  reduce  to. 


/  0u\  .  9u2  ^,x*+o 

'•"■•'ax’;  •  '“mW- 

dux  0u  >  xp+0 


()u 


C33 


3  I 4 -a 


~  0 

-  0 

=  -QS(xi)F(t) 


(12a) 


(12b) 


(12c) 
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Examination  of  the  relations  (12a,b)  reveals  that,  if  (c.,4c^s  —  djs )  <loes  not 
vanish  (  which  is  the  case  ),  then  u \  and  it,  mnst  be  at  least  constant  and. 
in  accordance  with  (10),  equal.  We  choose  this  constant  as  zero  and  tints 
conclude  that  u,  and  tq  vanish  at  the  interface.  Finally,  we  satisfy  the  normal 
stress  discontinuity  condition  (13c)  by  assigning 

C33~lx?+o  =  -QS{xx)F(t)/ 2,  c3 3“  =  Q6{x,)F(t)l 2  (1.13) 

<Xr3  1  c)x-i  ! 

Collecting  the  above  conditions  we  finally  summarize  the  condition  at  the 
artificial  interface  as 


it!  =  0 
u,  =  0 

tii  =  o 
u2  =  0 


for  ,r3  >  x3  at  x3  =  x3 


for  x3  <  x3  at  x3  =  xp 


(1.15) 


1.4  Formal  solution 

Following  the  procedure  of  [1],  we  outline  the  steps  leading  to  a  formal  solutions 
of  Eqs.(8a,b,c)  for  each  of  the  two  semi-spaces.  Since  the  body  force  has  been 
replaced  by  a  ’’boundary”  condition,  we  drop  /,  from  Eqs(8a.b,c).  We  then 
assume  harmonic  solutions  followed  by  applying  the  Fourier  transform  to  these 
equations  in  accordance  with 


r  - 

J  -oo 


-*r'dx  i 


(16a) 

(16b) 
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The  general  solution  of  the  resulting  differential  equations  is  then  sought  in 
the  form 


u,  =  r,e-a*'3,  1  =  1,2,3  (1.17) 

leading  to  the  characteristic  equation 


An 

A, 2 

a,3\ 

(o\\ 

A  ,2 

A  22 

A  23 

(’>  =  0 

(i 

A 13 

•^■23 

A  . 33  / 

\  l  3  / 

where  the  various  entries  A,:  are  given  in  Appendix  A.  Note  from  (18)  that 
the  A tJ  matrix  is  symmetric. 

For  the  existence  of  nontrivial  solutions  in  L't  the  determinant  in  Eq.  (18). 
must,  vanish  giving  an  algebraic  equation  relating  a  to  a;.  This  is  obviously  an 
alternative  presentation  of  Christoffe!  equation[l].  The  difference  is  that  we 
are  now  solving  for  a  in  terms  of  w  as  compared  with  solving  for  the  phase 
velocity  for  a  given  propagation  direction.  Setting  the  determinant  equal  to 
zero,  one  obtains  a  sixth  order  equation  in  a  (cubic  in  a2)  which  is  written 
symbolically  as 

A! a6  -1-  A2a4  +  A3q2  +  .4^  =  0  (1.191 

with  its  coefficients  are  given  in  Appendix  B.  Eq.(19)  admits  six  solutions  for 
a.  These  o's  have  the  further  properties  that 

ft 2  =  —a,,  a4  =  o6  =  -«5  (1-20) 


Furthermore  for  each  or,  equation!  18)  yields  the  displacement  amplitude  ratios 
Vy  —  L  2<j/f '  1  q  i  tWg  L  2q /  F  1 q 


A  1  I  A  23  —  A  12  A  13 

A 12  A  23  ~  A  22  A 13 
A  j  2  A  23  —  A  13  A  2  2 

A23A23  ~  A22  A;s;j 


(21a; 


(21b) 


II 


Finally,  using  superposition,  we  write  the  formal  solutions  for  the  displace¬ 
ments  of  equations  (8a,b,c)  and  their  associated  stress  components  using  Kq.(6) 


(«i,  u2,  u3)  =  £(1-'V  «•’,)£ :i,e  -c?)  ( 1 .22 ) 

t?=l 

(0-33^13.^3)  =  't(D]vD2g.D3q)Uue-*^-IV  (1.23) 

?=t 

where 

Dig  =  j£{Ci3  +  C360,)  ~  CxiCtgWg  (1.24) 

Dig  =  c$5(jtWg  -  Og)  -  cA5aqVg  ( 1 .25) 

Dig  =  CgsU^Wg  -  aq)  -  CggOtgVq,  </  =  1 .2,. .,6  (1.26) 


With  reference  to  equation  the  relation  (2d)  and  to  the  \v  entries  of  Appendix 
A  and  by  inspection  of  equations  (20)  and  (24-26).  one  recognizes  the  relations 
(see  Ref. [2]) 


V2 

=  vu 

u4 

-  V3, 

ve 

=  y.s 

U?2 

=  -Wi, 

Wg 

II 

1 

§ 

w6 

=  -lt?5 

D\2 

—  £11, 

D\4 

—  D 13, 

D\e 

=  D 15 

D22 

—  —D  21, 

D24 

—  —D23, 

D26 

=  —  £^25 

D32 

—  ~D3 1, 

D34 

—  —  Du, 

Die 

=  Die 

The  above  solutions  with  their  various  properties  can  now  be  specialized  to 
both  artificial  half-spaces  by  the  following  steps.  Inspection  of  the  above  so¬ 
lutions  indicate  that  each  consists  of  three  pairs  of  wave  components,  each 
pair  propagating  in  mirror  image  fashion  with  respect  to  the  interface,  namely 
along  positive  and  negative  13-directions.  Since  propagation  is  expected  to 
emanate  from  the  interface  into  both  media,  we  arbitrary  reserve  <71. 73,  and 
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f/5  for  the  lower  half-space;  the  remaining  one's,  namely  described  with  (/_>.</ 
and  (j( i  for  t he  upper  one.  We  list  the  formal  solution  in  the  lower  and  lower 
half-spaces  according  to 

(<"'1,  u  > ,  u-s)  =  Y  t’,  -  )/  igf- 

ij=l.3.S 

(^33-^13.-5-23)  =  Y 

q—  l  ,3,5 

(til.  1*2,63)  =  Y  (1.  i’,,te?)t;i,c-,*,<X;,-'r^) 

<}=2.4.t> 

( 5.33- 5,3.  <r23)=  Y 

<7=2,4.6 

1.5  Solutions  in  infinite  media 

At  this  point,  we  have  presented  a  formal  solution  of  the  field  equation  in  a 
monoclinic  medium.  The  amplitudes  L\n  are  the  unknowns.  The  amplitudes 
b\g  will  now  be  determined  by  implementing  the  artificial  interface  conditions 
(14-15).  To  this  end,  if  (28)  is  subjected  to  the  conditions  (14)  and  (20)  to 
(15),  one  finally  solves  the  displacement  amplitudes  as 

C:u  =  —Uu  —  (u.5  —  v3)Q /{'2cj3Dum)  (30a) 

b\z=—U\4  —  (t’i  —  v*,)Q  /  ('IcyzD  um )  (30b) 

i  1 5  =  U\ 6  =  (t'3  —  vi)Q/{2cxiDum)  (30c) 

where 

Dum  =  vi{a3w3  -  atsWs)  +  v3{a5w3  -  a,te,)  -f  !;5(a,u;i  -  a3w3)( !  .31) 

It  is  interesting  to  note  that  Dum  =  0  defines  an  equivalent  Christoffel  charac¬ 
teristic  equation  for  the  propagation  of  bulk  waves  in  the  medium  (-].  With 


* 


these  solutions  for  the  wave  amplitudes,  solutions  in  the  upper  region  can  be 
written  in  terms  of  q  —  1,3,5  as 

(«i,  »2,  un)  =  (-1,  -i’.P  wq)Uiqe~a‘,ljr'~jr') 

(*33, *13, *23)=  E  (-Dlq.D.2q,D3q)l\,e-^-^  i!  '52t 

7=  1,3,5 

In  summary,  solutions  (28)  and  (32)  with  the  amplitude  solutions  (30)  uniquely 
define  the  propagation  fields  in  the  lower  and  upper  "artificial"  semi-spaces. 
In  other  words,  their  combination  constitute  the  total  solutions  for  the  infinite 
medium. 


1.6  Solution  in  half-spaces 

We  now  adapt  the  solutions  of  the  infinite  media  (28)  and  (32)  to  solve  for 
the  case  where  the  free  boundary  is  intercepts  the  propagating  pulse  at  some 
arbitrary  location  parallel  to  the  plane  x3  —  0.  We  assume  that  the  free 
boundary  is  located  at  x3  =  —d  as  depicted  in  figure  2.  This  implies  that 
the  free  boundary  is  located  in  the  upper  region  and  thus  can  only  interfere 
with  the  propagation  fields  in  the  negative  x3  direction.  For  this  case,  the 
solution  (32)  will  constitute  an  incident  wave  on  free  surface.  As  a  results 
waves  will  reflect  from  the  free  boundary  and  propagate  in  the  positive  x3 
direction.  Thus,  appropriate  formal  solutions  for  the  reflected  waves  can  be 
adapted  from  the  solution  (28)  in  accordance  with  (Note  nowr  that  x3  does  not 
appear  because  solutions  are  referred  to  the  origin  r3  =  0) 


“21  “3)  =  E  (1  .vq,wq)Ulqe 

(1.33) 

17=1,3.5 

(1-34) 

17=1,3.5 
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With  this  the  total  solution  for  the  semi-spare  (designator!  with  superscript 
s')  which  is  required  to  satisfy  the  stress  free  boundary  condition  is  obtained 
by  superposing  the  incident  waves  and  reflected  waves  in  accordance  with 


+  if  (i.'V 

<j=  1 .3.5 

+  Z\D[,.D2,.D„)U;,e~"'' 

9=  1.3.5 

The  boundary  condition  is  given  by 


(1.35) 


(1.36) 


~S  ~  5  -5 

°I3  a23  —  a33 


0,  at  x3  =  —d 


By  imposing  the  boundary  conditions  ( 37 )  on  Eq.(6).  linear  system  is  obt  ained 
by 


where 


D\\  DIZ  £)15 

D2\  D  23  D  25 

^31  ^33  D  35 


^,3^3 

^5 


(1.38) 


/?i  = 

/?3  = 

As  = 


-(f/ffDnET  +  +  U$DlsEl) 

(UflDnEl  +  t/gflofS  +  U$D25El) 

{USDzxE[  +  i/g#^  +  U^DkED 


(39a) 

(39b) 

(39c) 


Using  the  standard  Cramer’s  rule  we  solve  the  reflected  amplitudes  as 


DsmE{ 

-1 

D  „ri  E-\ 


(RiGu  —  R3G  2t  +  Rf>G^\ ) 

( R\  G 13  —  R3G23  +  R-c,  C33 ) 


(40a) 

(40b) 
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rr 

L  IS 


Dsm  E$ 


{R\G\5  —  R3G0 5  +  R'X j'  jsj 


(40c) 


with 

£,  =  ea*d,  E;  =  e-°^+dK  q- 1,3.5  (1.41) 

and 

Dsm  =  DnG\\  —  D21G21  +  DiiG'ji  (1.42) 

where  Gqr  are  given  in  Appendix  B.  It  is  worth  noting  that  D,m  =  0  defines 
the  characteristic  equations  for  the  propagation  of  Rayleigh  (surface)  wave  on 
the  free  surface  [1] 


1.7  Solution  for  infinite  plates 

We  now  consider  the  response  to  the  line  load  in  an  infinite  plates  having  the 
thickness  2d.  Since  the  elastic  waves  in  the  plates  are  reflected  at  both  free 
boundaries  x3  =  d  and  x3  =  —  d  as  shown  in  figure  3,  the  solutions  are  obvi¬ 
ously  more  complicated  than  for  the  half  space  case,  we  however,  the  solution 
procedures  are  the  same.  To  this  end,  we  adapt  the  solutions  of  the  infinite 
space  (28,32)  in  order  to  solve  for  the  case  where  the  both  free  boundaries  in¬ 
tercept  the  propagating  field.  This  implies  that  the  free  boundaries  are  in  the 
upper  and  lower  semi-space  and  interfere  with  the  propagation  in  the  negative 
and  positive  x3  directions.  For  this  case,  the  both  upgoing  and  downgoing 
waves  in  the  infinite  space  will  constitute  incident  waves.  As  a  result  upgoing 
waves  will  reflect  from  the  free  boundary  x3  =  —  d  and  propagate  in  the  posi¬ 
tive  x3  direction  whereas  downgoing  waves  will  reflect  from  the  free  boundary 
x3  =  d  and  propagate  in  the  negative  x3  direction.  Thus,  total  formal  solu¬ 
tions  for  the  plate  (designated  with  superscript  ’p‘)  can  be  adapted  from  the 
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solution  for  the  infinite  space  and  free  waves  (i.e.  the  six  scattered  ones)  in 
accordance  with 


for  the  upper  artificial  layer  ( —d  <  x3  <  X3) 

(a?,as,s5)  =  E 

1?=  1.3,5 

9  =  1 

9=  1,3.5 

9=1 

for  the  lower  artificial  layer  (x^  <  X3  <  d) 

9=  1.3, 5 

+E(i,»„ul,)c;,e—« 

9=1 

9=1,3, 5 

+E(£>i„D2„Z?3,)t/1r,c-a*xs 

9=1 


(1.43) 


(1.14) 


(1.45) 


(1.46) 


The  polarized  amplitude,  f/,r  ,  is  again  determined  by  applying  the  free  traction 
on  the  free  surfaces  given  by 

^33  =  £13  =  ^23  =  0,  at  x3  =  id  (1.47) 

Now,  substituting  the  solution  for  the  plate  (??,??)  into  (??)  with  the  aid  of 
(27)  yields  a  linear  system  to  be 


D\\E\ 

D13E3 

DisEs 

Du  Ei 

D13E3 

Di$Es 

d21ex 

D23E3 

D25  E$ 

—D21E1 

D23E3 

-D25Ej 

D31  E\ 

D33E3 

D35  £5 

—D31  Ei 

D33E3 

—  D35  Ei 

DUE 1 

D13E3 

D15E5 

Du  Ei 

D13E3 

D\ 5  E$ 

-D2iEj 

—  D23E3 

—  D23E5 

D21E1 

D23E3 

DisEs 

-D-i\  Ei 

—  D33  E3 

—  O35  £V, 

D31  Ei 

D33E3 

D35E5 

(1,48) 
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where,  besides  the  definitions  of  R j,  R3,  and  R$  of  (39). 


r2  = 

(U"DUE\  +  u$dX3e;  +  u$d15e:) 

(49a) 

ra  — 

(u[\]d2Xe\  +  u$d23e;  +  l’SdkED 

(49b) 

Re  = 

(U$D3}E\  +  U$D33Et3  +  uttDxEl) 

(49c) 

with 

K  = 

e~aqd,  Er  =  e-a'^rf-x3l,  q=l,3,5 

(1.50) 

After  straight  forward  algebraic  manipulation,  the  amplitude  U\q  is  algebraically 
simplified  to 


UTU  =  -UTU 
3  =  -uu 

U15  =  -Ul 6 


<? 


sinh(«3d)  sinh(o'5(i) 


[(/2n^ll  —  /?22#12)(#52^41  ~ 


Dpmcoth(a5d)£)f5 

—  (#33#41  ~  ^?44  R-i2){  R^22  R\\  ~  ^I2^2l)} 
^sinh(a,</)sinh(a5d)r,0  „  D  „  u  u  u 

-V-7; - . ,xn2  1^11^21  -  ri22n22){nA2n5i 

Dpm  coth (a5d)Df5 

“*** { R33R51  -  R«Hs2){HnH2X  -  Hu  H22)] 

sinh(aid)  sinh(a3<f) .  p  n  uu  u 

— 7Tn2  llb-u  H31  -  R-22n3x){tiA2ne\  ~ 
Dpn  coth(a3a)£>f3 

—  (^33^61  —  f?44  #62X^12^31  ~  #11^32)] 


//51//42) 

(51a) 

(51b) 

R  41  ^62) 

(51c) 


with 

£)pm  =  — 4  cosh(aid)  cosh(a3d)  cosh(Qsd)x 

[Z?nGu  coth(ajd)  -  Di3Gt3coth(a3<f)  +  AsGis  coth(a5d)]x 
[DnGu  tanlXcqd)  -  DX3GX3  tanh(a3<f)  -f  DXSGX5  tanh(a5d)]  ( 1 .52) 

wher^  Rkk  and  H^i  are  given  in  Appendix  B.  Note  that  Dprn  =  0  defines  the 
characteristic  equations  for  the  propagation  of  free  wave  in  the  plate. 
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1.8 


Numerical  results  and  discussions 


In  this  section  we  present  numerical  illustrations  of  the  above  analysis.  We 
choose  for  our  illustration  a  graphite-epoxy  composite  material.  The  mi¬ 
crostructure  of  this  particular  material  has  been  extensively  studied  and  mod¬ 
els  predicting  its  effective  anisotropic  properties  are  available  elsewhere  [l]. 
Based  upon  a  60%  graphite  fiber  volume  fraction  these  effective  properties  are 
given  with  respect  to  the  reference  coordinate  system  x'  as 


/  155.60  3.70 

3.70  15.95 
3.70  4.33 

0  0 

0  0 

\  o  o 


3.70  0  0 

4.33  0  0 

15.95  0  0 

0  5.81  0 

0  0  7.46 

0  0  0 


0  \ 
0 
0 
0 
0 


x  1010 


7.46  / 


[  N/m2 


with  p  =  1.6  g/cm3.  For  a  rotation  of  <f>  =  45°,  for  example,  these  properties 
transform  to 


/  51.20 

36.28 

4.03 

0 

0 

36.13  \ 

36.28 

51.20 

4.03 

0 

0 

36.13 

4.03 

4.03 

16.00 

0 

0 

0.31 

0 

0 

0 

6.64 

0 

0 

0 

0 

0 

0 

6.64 

0 

V  36.13 

36.13 

0.31 

0 

0 

40.02  j 

which  confirms  the  earlier  conclusion  that  the  transformed  matrix  takes  the 
format  of  monoclinic  symmetry.  Having  chosen  the  material  we  now  sum¬ 
marize  a  Tow  chart”  like  procedure  for  our  subsequent  calculations.  After 
specifying  the  azimuthal  angle  <j>  (namely,  the  line  load  direction)  we  proceed 
in  the  first  step  to  evaluate  a' s  for  given  values  of  the  wavenumber  if.  We  then 
sort  out  the  various  a’s  in  accordance  with  the  required  format  of  equation 
(28-29).  In  the  second  step  we  evaluate  the  displacement  ratios.  By  now,  we 
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are  ready  to  calculate  the  inverse  transforms  in  accordance  with 


1  t°° 

«,  =  —  /  u<ejir'dZ,  i  —  1,2.3  (1.53) 

Depending  upon  whether  u,  is  symmetric  or  antisymmetric,  computational 
effort  would  be  greatly  reduced.  We  first  would  like  to  present  solutions  to 
symmetry  situations.  If  aq  is  symmetric  with  respect  to  £.  Then  vq  will  be 
symmetric  while  wq,  Diq ,  Dsm,  and  Dpm  will  be  antisymmetric  with  respect 
to  <f.  It  then  follows  that  up  and  u2  are  symmetric  while  u3  is  antisymmetric. 
Accordingly  equation  (??)  simplifies  to 
1  r°° 

u,  =  -  ;u,sin(£xi)d£,  (  =  1.2  (54a) 

7T  J  0 

1  f°° 

u3  =  -  /  U3Cos(^!)d^  (54b) 

7T  Jo 


The  calculations  of  displacement  and  stress  components  are  then  formed  straight¬ 
forward,  except  at  singularities.  These  singularities  are  poles  corresponding 
to  the  zeros  of  Dsm  and  Dpm,  namely,  the  characteristic  equations  for  a  sur¬ 
face  wave  in  the  halfspace  and  free  waves  in  the  plate,  respectively.  We  find 
that  the  poles  do  not  exist  in  the  second  and  the  fourth  quadrant  of  figure  4, 
namely  for  Re(£)  >  0  and  /m(<f)  <  0  and  that  Dum  of  the  infinite  medium 
has  no  poles.  Thus,  u,  in  the  infinite  space  is  integrated  along  the  real  £-axis. 
Two  different  procedures  have  been  used  in  dealing  with  these  poles  in  the 
integration  process.  The  first  requires  the  removal  of  these  singularities  from 
the  integrals.  The  second,  based  upon  the  work  of  Neerhoff  et  al  [18].  merely 
deforms  the  contour  of  integration  below  the  real  (f-axis,  as  shown  in  figure 
4,  so  that  no  poles  occur  on  the  path  of  integration.  In  order  to  optimize 
accuracy,  judicious  choices  of  the  parameter,  are  required.  The  value  of 
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imaginary  part  of  needs  to  be  large  enough  to  offset  the  contributions 
of  singular  points,  but  not  two  large  since  it  reduces  the  accuracy.  In  our 
calculations,  we  choose  the  value  <f,  =0.1.  For  our  numerical  illustrations 
we  arbitrary  choose  the  azimuthal  angle  0  =  45°  and  the  angular  frequency. 
u>  —  3.0  x  106  rad/sec. 

Figures  5-7  depict  the  displacement  and  a  normal  stress  components 
Ui,u2> ^3, <?33  in  the  media  at  the  location  x3  =  —3  mm.  The  location  of 
the  line  load  is  chosen  at  the  origin  of  the  x ,  coordinate  system  for  all  three 
geometries  (infinite-  and  half-space,  plate):  thus  ut  and  <733  are  symmetric 
with  respect  to  x3  axis  while  Uj  and  u2  are  antisymmetric.  Note  from  the  fig¬ 
ures  that  the  amplitudes  of  the  displacement  and  stress  component  are  slowly 
decreasing  with  oscillations  for  a  long  distances  from  the  location  of  line  load 
in  the  infinite  and  semi-infinite  space.  However,  the  amplitudes  in  the  plate 
are  oscillating  along  the  location  x3  =  —3  mm  with  persisting  amplitudes. 

Figures  8-10  present  the  displacement  mode  (u],u3)  at  various  locations. 
We  calculate  deformation  fields  at  the  locations  designated  by  broken  various 
radii  (2.4mm,  5.55mm,  and  7.2mm  in  infinite  media;  2.4mm  and  5.55mm  in 
infinite  media;  and  4.0mm  in  the  plate)  away  from  the  source  location.  For 
these  investigations  the  line  load  is  located  at  the  origin  of  x,  coordinate  system 
in  the  infinite  and  half-space  media,  and  at  (0,0,1)  in  the  plate.  Solid  curves 
represent  deformations  defined  by  xj  +  uj  x  102,  x3  +  u3  x  102.  Note  that  in  the 
infinite  media,  u\  is  zero  on  the  xt  axis  and  that  u3  has  maximum  amplitude 
on  the  X3.  However,  in  the  semi-infinite  space  and  the  plate  the  displacement 
fields  are  much  affected  by  the  reflected  waves.  The  displacement  amplitude 
in  the  semi-infinite  space  is  seen  to  be  greater  than  that  in  the  infinite  space. 
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Figure  1.1:  A  applied  line  load  in  an  anisotropic  infinite  media 


Figure  1.2:  Semi-infinite  media 
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Figure  1.3:  The  infinite  plate 
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Figure  1.4:  Complex  integral  contour 
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Figure  1.5:  Displacements  and  stre: 
x3  =  -3  mm  :  (a)  (b)  u2,  (c)  w3, 


UjXlO2  UjXlO 


-9.0  -6.0  -3.0  0.0_  3.0  6.0  9.0 

X„  (X1+U,X102) 


Figure  1.8:  Displacement  modes  in  the  infinite  space  of  graphite-epoxy.  Bro¬ 
ken  circles  are  chosen  as  receiver  location  (xl7x3)  and  solid  curves  are  displace¬ 
ments  modes  (x,  +  U,  x  102,x3  4-  u3  x  102)  with  respect  to  broken  circles.  The 
line  load  is  located  at  (0,  0,  1) 
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-6.0  -3.0  0.0  3.0  6.0 

X,,  (x,+u,x102) 


Figure  1.9:  Displacement  modes  in  the  half-space  of  graphite-epoxy.  Broken 
circles  are  chosen  as  receiver  location  (x],x3)  and  solid  curves  are  displace¬ 
ments  modes  (xt  +  u\  x  102,x3  -j-u3  x  102)  with  respect  to  broken  circles.  Free 
surface  is  x3  =  —5.55  mm  and  the  line  load  is  located  at  the  origin  of  the 
system. 
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-6.0  -4.0  -2.0  0.0  2.0  4.0  6.0 

x,,  (x,+u,x10z) 


Figure  1.10:  Displacement  modes  in  the  plate  of  graphite-epoxy.  Broken  circle 
is  chosen  as  receiver  location  (xl5x3)  and  solid  curve  is  displacement  mode 
(xi  +  U\  x  102,x3  +  u3  x  102)  with  respect  to  broken  circle.  Free  surface  is 
x3  =  ±5.55  mm  and  the  line  load  is  located  at  (0,  0,  1). 
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Appendix  A 

Various  coefficients  of  characteristic  eq.(14)  are  given  by 
An  =  e55Q2  —  Cn£2  +  puj2 

A 12  =  C45«2  —  Ci6<£2 

Ai3  =  -j£o(ci3  +  c55)  (A.l) 

A22  =  C44O!2  -  C66£2  +  pu>2 

A  23  =  -j£o(C; 36  +  C45) 

A33  =  C33«2  —  0551^ 2  +  pOJ2 


At  —  C33C25  —  C33C44C55 
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A-2  =  [C11C33C44  +  C55  F\  —  C16f33C45  —  C45  F>  +  (Cj3  -f  C55  )  Ojjif  2~r 

[c-55  ( C33  +  C44)  +  C33C44  —  C'i5\pU! 

A3  =  [(Ci3  +  C55)F4  +  C45C16C55  +  C16F2  —  C^Ogg  —  Cn 

[(c13  +  c55)2  +  ~c45c16  —  F\  ~  055(055  +  Cge)  —  Ct  1(033  +  C44)]pu,’£2  — 

(033  +  C44  +  C55)p2u;2  (A. 2) 

A4  —  (O11O55O66  Ojg055)^  “f“  [011(055  4"  Ogg  A 

(Oil  +  C55  +  C6e)p2u>2£2  +  p3u>3 

F\  —  C33C66  +  C44C55  —  (C36  +  C45)2 

F 2  —  c33c16  +  c45c55  ~  (c36  +  c45)(ci3  +  C55) 

O3  =  045(036  +  C45)  —  044(013  +  C55)  (A. 3) 

F 4  =  C66(Ci3  +  C55)  —  Oi6(c36  +  C45) 

Appendix  B: 

Various  coefficients  of  eq. (32-40)  are  given  by 

(?n  =  O23O35  —  0,33/225 

G21  =  O13O35  —  D33O15 

G31  =  O13O25  —  O15O23 

f?13  =:  D31D25  —  O21O35 

G23  =  O31O15  —  O11O35  (B.l) 

G33  =  O15O21  —  O11O25 

G'15  =  O21O33  —  O31 O23 

G25  =  On  O33  —  O31 013 

G35  =  O11O23  — O13O21 

//n  =  O15O23 coth(a,5(f)  -  O13O25 coth(o3<i) 


3.3 


Hu  =  #is#33  coth(asd)  —  #i3#35  coth(a3<i) 

# 21  =  DxsD2i  coth(a5d)  -  #u#25  coth(axd) 

#22  =  #15  D31  coth(a.'5</)  —  DxlDMcoth(ald) 

#31  =  #i3#2i  coth(of3<f)  —  D\\ Z?23  coth(oti(Z) 

#32  =  #i3#3i  coth(c*3d)  —  Dn #33 coth(o,i d) 

#41  —  D15D23  coth(a3tf)  —  D13D25  coth(o,5d)  (B.2) 

#42  =  £>15^33  coth(a3^)  -  #i3#35  coth(a5d) 

#51  =  #is#2i  coth(ai^)  -  #n#25  coth(a5.f) 

#52  =  Di5#3i  coth(ai<i)  —  #11  #35  coth(o,5<f} 

#61  =  #i3#2i  coth(airf)  -  DnD23coth(a3ii) 

#62  =  #i3#3i  coth(airf)  -  #u #33  coth(a3d) 

and 

#11  =  (R3  -  R6)Dlscoth(a5d)  -  (Rx  +  R4)D35 

#22  =  (#2  -  #5)  #15  coth(o5d)  -  (#1  4*  #4)  #25 

#33  =  (#1  —  #4)#35  c°th(a5<f)  —  (#3  4-  /?e)#15 

#44  =  (#1  —  #4)#25  coth(osd)  —  (#2  +  #5)#15 

#55  =  (#3-«6)#l3COth(a3d)-(#l+#4)#33  (B.3) 

#66  =  (#2  ~  #5)  #13  COth(Q3d)  —  (#1  +  #4  )#23 

#77  =  (#1  —  #4)#33  COth(a3d)  —  (#3  +  #6)#13 

#68  =  (#1  —  #4) #23  COth(Of3(/)  —  (#2  +  #5) #13 
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Chapter  2 

WAVE  PROPAGATION  IN 
ANISOTROPIC  MEDIA  DUE 
TO  TRANSIENT  SOURCES 


2.1  Introduction 

Plane  harmonic  wave  interaction  with  homogeneous  elastic  anisotropic  media, 
in  general,  and  with  layered  anisotropic  media,  in  particular,  have  been  ex¬ 
tensively  investigated  in  the  past  decade  or  so.  This  advancement  has  been 
prompted,  at  least  from  a  mechanics  point  of  view,  by  the  increased  use  of 
advanced  composite  materials  in  many  structural  applications.  Being  both 
anisotropic  and  dispersive,  composite  materials  required  in  depth  understand¬ 
ing  of  their  mechanical  behavior.  The  list  of  relevant  literature  is  rather  long 
and  we  thus  refer  the  reader  to  selective  recent  works  for  further  references 
[1-6]. 

The  difficulties  inherent  in  the  treatment  of  wave  propagation  in  anisotropic 
media  in  general  can  be  illustrated  by  considering  slowness  wave  surface  meth¬ 
ods  for  infinite  media  [1,7-9].  In  this  analysis,  the  slowness  surface  of  an 
isotropic  material  consists  of  two  concentric  spherical  sheets,  an  inner  one  rep- 
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resenting  the  longitudinal  wave  mode  and  an  outer  one  representing  the  two 
degenerate  transverse  modes.  In  the  anisotropic  case,  we  find  three  general 
wave  surfaces,  one  for  the  so-called  quasi-longitudinal  wave  and  two  nonde¬ 
generate  sheets  for  the  quasi-transverse  waves.  Moreover,  the  surfaces  will  no 
longer  be  spherical  in  shape  but  will  reflect  the  elastic  symmetry  or  asymmetry 
of  the  material. 

Perhaps  the  most  severe  consequence  of  elastic  anisotropy  in  infinite  media 
is  the  loss  of  pure  wave  modes  for  general  propagation  directions.  This  fact 
also  implies  that  the  direction  of  energy  flow  (i.e.,  group  velocity)  does  not 
generally  coincide  with  the  wave  vector,  or  wavefront  normal.  Thus,  uncoupled 
pure  potentials  (such  as  are  found  in  the  isotropic  case)  are  much  easier  to 
treat  than  the  mixed  modes  characteristic  of  anisotropic  materials.  For  wave 
propagation  in  directions  of  symmetry  some  wave  types  revert  to  pure  modes, 
leading  to  a  simpler  characteristic  equation  of  lower  order. 

A  key  condition  which  was  found  to  facilitate  the  previous  analysis  is  that 
at  any  boundary  all  wave  vectors  must  lie  in  the  same  plane.  This  requirement 
implies  that  the  response  of  the  media  will  be  independent  of  the  in-plane 
coordinate  normal  to  the  propagation  direction.  Accordingly  the  analysis  was 
conducted  in  a  coordinate  system  formed  by  the  line  load  direction  and  its 
normals  rather  than  one  determined  by  material  symmetry  axes.  This  choice 
leads  to  a  significant  simplification  in  algebraic  analysis  and  computations. 

In  comparision  to  the  extensive  literature  on  the  interaction  of  plane 
harmonic  waves  with  anisotropic  media,  very  little  work  is  available  on  the 
response  of  such  media  to  concentrated  source  loading.  Here  concentrated 
sources  include  point  as  well  as  line  loads  of  which  harmonically  pulsating  and 
transient  sources  are  common  types  of  such  loading.  Understanding  the  re¬ 
sponse  of  elastic  solids  to  internal  mechanical  sources  has  long  been  of  interest 


to  researchers  in  classical  fields  such  as  acoustics,  seismology,  as  well  as  modern 
fields  of  application  like  ultrasonics  and  acoustic  emission.  It  is  known  that 
whenever  a  material  undergoes  a  local  failure,  elastic  waves  are  generated  due 
to  the  rapid  release  of  localized  strain  energy.  Such  radiation,  for  example,  is 
known  as  acoustic  emission  in  the  field  of  nondestructive  testing  of  materials. 
In  seismology  it  is  of  course  known  as  earthquake. 

A  quick  review  of  available  literature  on  this  subject  reveals  that  most  of 
the  work  done  so  far  is  carried  out  on  isotropic  media.  The  effect  of  imposed 
line  load  in  homogeneous  isotropic  media  has  been  discussed  by  several  in¬ 
vestigators  ever  since  Lord  Rayleigh  discovered  the  existence  of  surface  waves 
on  the  surfaces  of  solids  [10].  An  account  of  the  literature  dealing  with  this 
problem  through  1957  can  be  found  in  Ewing,  Jardetsky  and  Press  [18].  Most 
of  the  earlier  work  [13-15]  followed  Lamb  [11,12],  who  apparently  was  the  first 
to  consider  the  motion  of  semi-infinite  space  caused  by  a  vertically  applied 
line  load  on  the  free  surface  or  within  the  medium.  It  was  shown  that  dis¬ 
placements  at  large  distance  consists  of  a  series  of  events  which  corresponds  to 
the  arrival  of  longitudinal,  shear,  and  Rayleigh  surface  wave.  The  analytical 
approach  used  in  the  above  mentioned  investigations  and  others  [13-15]  can 
be  summarized  as  follows  :  The  steady  state  problem  for  harmonically  pul¬ 
sating  source  in  infinite  isotropic  media  is  solved  at  first  and  then  generalized 
to  the  case  of  half-space  using  superposition  technique.  For  transient  source 
loading  results  can  be  obtained  from  those  corresponding  to  harmonic  ones  by 
a  Fourier  integral  approach.  The  resulting  double  integral  could  be  evaluated 
only  by  considering  large  distances. 

A  modern  alternative  approach  has  been  suggested  by  Cagniard  [16].  He 
showed  that  a  suitable  deformation  of  the  integral  contour  not  only  resulted  in 
considerable  analytical  simplification  but  led  to  exact,  closed  form,  algebraic 
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expressions  for  the  displacements  as  functions  of  time.  Subsequently,  many 
investigators  have  obtained  the  disturbances  due  to  sources  in  the  isotropic 
media  by  applying  Cargniard  approach  [17].  DeHoop  simplified  Cargniard's 
method  considerably  and  applied  to  cases  involving  impulsive  line  and  point 
sources  in  infinite  homogeneous  isotropic  media  [18-21].  Kraut(196d)  ap¬ 
plied  the  Cargniard-DeHoop  method  for  investigating  the  transient  distur¬ 
bance  caused  by  a  surface  line  load  in  the  anisotiopic  half-space.  Van  der 
Hijden(1987)  investigated  the  features  of  the  wave  propagation  in  infinite 
anisotropic  media  generated  by  a  mechanical  line  source  [28]. 

In  this  paper,  we  closely  follow  the  formal  developments  in  previous  works 
[i-6]  and  study  the  response  of  two  anisotropic  systems  to  transient  buried  line 
loads.  This  analysis  include  infinite  and  semi-infinite  systems.  The  foregoing 
illustration  shall  be  carried  out  on  ani  V.opic  media  possessing  monoclinic 
or  higher  symmetry.  The  load  will  be  in  the  form  of  a  normal  stress  load 
acting  at  an  arbitrary  direction  in  the  plane  of  symmetry  of  the  material.  One 
then  uses  a  building  block  approach  in  which  one  starts  by  driving  results 
for  an  infinite  media.  Subsequently  one  obtains  the  results  for  the  half-space 
employing  superposition  of  the  infinite  medium  solution  together  a  scattered 
soluti-  from  the  boundary.  The  sum  of  both  solutions  has  to  satisfy  stress 
free  boundary  conditions  thereby  yielding  to  complete  solutions.  Consequently 
explicit  solutions  for  the  particle  displacements  in  both  systems  are  obtained 
using  Cargniard-DeHoop  contour. 

This  work  will  accomplished  by  using  the  linear  transformation  approach 
in  which  one  identifies  the  line  lead  with  the  x2  direction.  This  implies  that,  all 
involved  field  variables  will  be  independent  of  the  x2  direction.  Nevertheless, 
and  in  general,  one  has  three  nonvanishing  particle  displacements.  Material 
systems  of  higher  symmetry,  such  as  orthotropic,  transversely  isotropic,  cubic. 
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Figure  2.1:  A  applied  line  load  in  an  anisotropic  infinite  media 

and  isotropic  are  contained  implicitly  in  the  analysis.  Numerical  results  are 
demonstrated,  drawn  from  concrete  examples  of  materials  belonging  to  sev¬ 
eral  of  these  symmetry  groups.  It  is  found  that  for  orthotropic  and  higher 
symmetry  materials  where  the  remaining  two  principal  axes  lie  in  the  plane  of 
symmetry,  the  particle  motions  in  the  sagittal  and  the  normal  to  it  uncouple 
if  propagation  occurs  along  either  of  these  in-plane  axes. 

2.2  Theoretical  development 

Consider  an  infinite  anisotropic  elastic  medium  possessing  monoclinic  symme¬ 
try.  The  medium  is  oriented  with  respect  to  the  reference  cartesian  coordinate 
system  x\  —  ,  x'2,  X3)  such  that  the  z’3  is  assumed  normal  to  its  plane  of 

symmetry  as  shown  in  figure  1.  The  plane  of  symmetry  defining  the  mono- 
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clinic  symmetry  is  thus  coincident  with  the  s\  —  s',  plane.  With  respect  to  this 
primed  coordinate  system,  the  equations  of  motion  in  the  medium  are  given 
by  [l] 


- 7  +  /,  =  P - 


/•) 


dx‘ 


dt> 


and,  from  the  general  constitutive  relations  for  anisotropic  media. 
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by  the  specialized  expanded  matrix  form  to  monoclinic  media 
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where  the  standard  contracted  subscript  notations  1  — >  11,  2  — *  22,  3  — *  33. 
4  — >  23.  5  — ! *  13,  and  6  -»  12,  to  replace  c^^i.j.kJ  =  1,2,3)  with  cvq[p.q  — 
1, 2, 6)  are  employed.  Here  cr'tJ ,  e'tj  and  u't  are  the  components  of  stress,  strain 
and  displacement,  respectively,  and  p'  is  the  material  density.  In  Eq.(1.3). 
7 ij  =  2e,j  (with  i  /  j  )  define  the  engineering  shear  strain  components. 

In  what  follows,  one  considers  the  response  of  the  infinite  medium  to  a 
uniform  transient  (time  dependent)  line  load  applied  along  a  direction  that 
makes  an  arbitrary  azimuthal  angle  0  with  the  x',  axis.  That  is,  the  direction 
0  =  0°  coincides  with  the  reference  coordinate  x\.  Since,  as  was  pointed  out 
in  the  introduction,  the  response  of  the  medium  to  such  a  wave  is  independent 
of  the  applied  line  direction.  The  analysis  is  conducted  in  a  transformed 
coordinate  system  x,  =  (xi,X2,x3)  formed  by  a  rotation  of  the  plane  x)  —  x'2 
through  the  angle  0  about  the  x!,  direction.  For  convenience  the  direction  x^ 
is  chosen  to  coincide  with  the  line  load  direction. 
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Since  CijU  is  a  fourth  order  tensor,  then  for  any  orthogonal  transformation 
of  the  primed  to  the  non- primed  coordinates  it  transforms  according  to 


C.jkl 


mnop 


(2.4) 


where  3t}  is  the  cosine  of  the  angle  between  x[  and  xt.  For  a  rotat  ion  of  ang! 
<?  in  the  x)  —  x'2  plane,  the  transformation  tensor  3,,  reduces  to 
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If  the  transformation  (1.5)  is  applied  to  Eq.(i.3).  one  obtains 
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where  the  relations  between  the  cvq  and  c'pq  entries  are  listed  in  [5 


>.6) 


Notice 
•  (1.3) 


that,  no  matter  what  rotational  angle  <}>  is  used,  the  zero  entries  in 
will  remain  zero  in  Eq.(1.6).  In  terms  of  the  rotated  coordinate  system  x*.. 
momentum  equations  can  be  written  in  the  form 


dtJij  d2ut 

'd^  +  f'  =  p-W 


(2.7) 


As  mentioned  earlier,  in  the  rotated  system,  the  elastic  wave  equations,  for 
wave  propagating  in  the  Xj  —  x3  plane,  are  independent  of  x2  .  Nevertheless, 
the  particle  motion  can  generally  have  three  nonzero  components  uj,u2,  and 
U3.  The  U2  displacement  can  be  identified  as  belonging  to  the  horizontally 
polarized  (SH)  wave.  Here,  the  equations  describing  these  three  wave  motions 
are  coupled  in  accordance  with 
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/,■  is  chosen  to  be  of  the  form 

/,•  =  Q,6(xj)S(x3  -  x§)F(t)  (2.9) 

2.3  Source  characterization 

Using  the  above  geometric  arrangement,  one  describes  the  dire  delta  line  load 
as  body  force  function  located  at  the  arbitrary  location  (x3  =  ).  This  choice 

of  source  location  is  arbitrary  and  is  not  necessary  for  the  cases  of  propaga¬ 
tion  in  either  the  infinite  nor  the  semi-infinite  media.  However,  as  will  be 
demonstrated,  this  choice  will  lead  to  great  simplification  in  the  algebraic  ma¬ 
nipulation  for  the  propagation  in  the  plate.  Solutions  to  the  present  problem 
can  be  accomplished  by  following  the  procedure  used  by  Achenbach  [32].  Ac¬ 
cording  to  this  procedure,  the  infinite  space  can  be  thought  of  as  consisting 
of  two  semi-spaces  whose  artificial  interface  contains  the  applied  load.  Thus 
the  upper  half-space  occupies  the  region  x3  >  x3  whereas  the  lower  half-space 
occupies  the  region  x3  <  x3.  As  a  result  of  this,  appropriate  condition  must 
be  specified  at  the  artificial  interface. 

In  order  to  be  able  to  specify  these  interface  conditions,  consider  a  very 
thin  ’’interface”  layer  extending  from  x3  =  x3  —  o  to  x3  =  x3  +  o  ;  Then 
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start  by  requiring  continuity  on  the  displacement  components  and  their  time 


derivatives,  namely 

u,(xi,X3.i)Q*0 

<?<n)Ui(Xi,X3,  t)  xp+0 

dtW  lr'-° 


0 

0,  i  =  1,2.3 


(10a) 

(10b) 


Integrating  Eqs.(8a.b,c)  across  the  interface  layer  and  using  the  continuity 
equations)  10a, b),  leads  to  the  discontinuity  in  the  displacement  spatial  deriva¬ 
tives  resulting  in  the  following  jump  conditions. 


dux 

{Css8x3 

■  r  <iU  :  '  I1!  —  J  _ 

- Qx6(Xl)F(t ) 

(11a) 

dux 

{c>sW3 

dll2  .Xj  +o 

-Q2S(Xx)F(t) 

(11b) 

du3  ,x^+o 

C33  77  rP_„  — 

dx3  0 

-QzHxx)F(t) 

(11c) 

Now,  consider  the  special  case  where  Q\  =  Q 2  =  0  and  Q3  =  Q.  For  this 


situation  equations  (11)  reduce  to. 

0u,  du2,,i?+o 

(CS5^  +  C45—  )  p  - 

dx3  dx3  x*  0 

dui  du2  xp+- 

(c,s— +  c„  — )U- 

du3 


c33 


1 3+0 


”3 

:JT?  +  0 


dx3 


<733|xp_ 


3TU  _ 


0 

0 

-Q8(xx)F(t) 


(12a) 


(12b) 

(12c) 


Examination  of  the  relations  (12a,b)  reveals  that,  if  (c44c55  —  C45)  does  not 
vanish  (  which  is  the  case  ),  then  u\  and  u2  must  be  at  least  constant  and,  in 
accordance  with  (10),  equal.  This  constant  is  then  chosen  to  be  zero  and  thus 
conclude  that  u \  and  u2  vanish  at  the  interface.  Finally,  satisfing  the  normal 

stress  discontinuity  condition  (13c)  and  assigning 

du  8u 

C33— =  -Q6(xl)F(t)/ 2,  C33— =  Q6(xl)F(t)l2  (2.13) 
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Collecting  the  above  conditions  one  can  summarize  the  condition  at  the  arti¬ 
ficial  interface  as 


«i  =  0 

U-2  =  0 

duj, 

UX-x 


=  --Q6(Xl)F(t) 


for  x3  >  x3  at  x3  =  x3  ('2.1 4  j 


«i  =  0 

14-2  —  0 
du3 

c33“^ 

OX-x 


-QS{x,)F(t) 


>  ,  for  x3  <  x3  at  x3  =  x3  (2.15} 


2.4  Integral  transforms  of  formal  solutions 


Following  the  procedure  of  [1],  let  us  outline  the  steps  leading  to  a  formal 
solutions  of  Eqs.(8a,b,c)  for  each  of  the  two  semi  spaces.  Since  the  body 
force  has  been  replaced  by  the  "artificial  interface-’  condition,  one  can  drop  /, 
from  Eqs(8a,b,c).  Assume  formal  solutions  followed  by  applying  the  Fourier 
transform  to  these  equations  in  accordance  with 


u,e  ptdt 
Ui€~}p7,Tl  dxi 


(2.16) 

(2.17) 


The  general  solution  of  the  resulting  differential  equations  is  then  sought  in 
the  form 


ut  =  Ute~pax\  i  —  1,2,3 
leading  to  the  characteristic  equation 

/  An  A12  A i3  \  /  U\  \ 

I  A 12  A22  A 23  1  I  U-i  I  =  0 

\  A 13  A 23  A  33/  \  U3  / 


(2.1S) 


(2.19) 


44 


where  the  various  entries  Au  are  given  in  Appendix  A.  Note  from  (1.19)  that 
the  A matrix  is  symmetric. 

For  the  existence  of  nontrivial  solutions  in  U,  the  determinant  in  Eq. 
(1.19),  must  vanish  giving  an  algebraic  equation  of  a.  This  is  an  alternative 
presentation  of  ChristofFel  equation  [2],  The  difference  is  that  one  now  solves 
for  a  as  compared  with  solving  for  the  phase  velocity  for  a  given  propagation 
direction.  Upon  setting  the  determinant  equal  to  zero,  one  obtains  a  sixth 
order  equation  in  a  (cubic  in  a2)  which  is  written  symbolically  as 

.4iq6  +  /Uo'4  -F  A.30^  +  A4  —  0  (2.20) 

with  its  coefficients  are  given  in  Appendix  B.  Eq.(1.20)  admits  six  solutions 
for  a.  These  a's  have  the  further  properties  that 

o2  =  -<*i,  =  -<*31  °6  =  -os  (2.21) 


Furthermore  for  each  a,  eq  jation(18)  yields  the  displacement  amplitude  ratios 

Vj  —  &2q / Wq  ~  ^3q/^-'\q 


4 1 1 A  23  —  A 12  A 13 

Ai2A23  —  A22Ai3 
A12A23  —  A13A22 

A23A23  —  A22A33 


(21a) 

(21b) 


Finally,  invoking  superposition,  one  can  write  the  formal  solutions  for  the 
displacements  of  equations  (8a,b,c)  and  their  associated  stress  components 
using  Eq.(  1.6)  as 

(ui,u2,u3)  =  £(l,u,,«t7)lV_pa,(a3-xJ)  (2.22) 

(d~33,  *3-13,^23)  =  y~!  P(Dui  Dig- fitalUlqe  pa,,(ri  (2.23) 

q= t 
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where 


D\q 

=  jl?(c,3  +  C36U,)  * 

-  c33aqwq 

(2.24) 

D2q 

=  c55(ji)ivq  -  aq) - 

-  c45  aqvq 

(2.25) 

Dzq 

=  c4b{jqwq  -  aq)  - 

~~  C^CkyVq  y 

sO 

II 

(2.26) 

With  reference  to  equation  the  relation  (20)  and  to  the  \t]  entries  of  Appendix 
A  and  by  inspection  of  equations  (20)  and  (24-26),  one  recognizes  the  relations 
(see  Ref. [2]) 


v2 

=  tq, 

va 

=  t’3, 

V6 

=  vb 

w2 

=  -Wi, 

Wq 

3 

1 

11 

we 

=  -iv5 

D\2 

—  Du, 

D\4 

=  D\ 3, 

D\e 

=  £>15 

d22 

=  —D2\, 

D2  4 

—  ~D2m 

£>26 

=  D2s 

D32 

—  —  An. 

Dm 

—  ~~d33, 

Dm 

—  —  D3S 

The  above  solutions  with  their  various  properties  can  now  be  specialized  to 
both  artificial  half-spaces  by  the  following  steps.  Inspection  of  the  above  so¬ 
lutions  indicate  that  each  consists  of  three  pairs  of  wave  components,  each 
pair  propagating  in  mirror  image  fashion  with  respect  to  the  interface,  namely 
along  positive  and  negative  redirections.  Since  propagation  is  expected  to 
emanate  from  the  interface  into  both  media,  we  arbitrary  reserve  q\,q3.  and 
q5  for  the  lower  half-space;  the  remaining  one’s,  namely  described  with  q2,q4< 
and  q6  for  the  upper  one.  We  list  the  formal  solution  in  the  lower  and  lower 
half-spaces  according  to 

(u1,u2,u3)=  Z  (l,vq,wq)Ulqe-paq{X3-x3] 

9=  1.3,5 

(<753,  V\Z,  <723)  =  Z  P(D\q,  Dlq,  D3q)i’iqe  po'^X3 

9=1,3, 5 

(ui,u2  u3)  =  Z  (1  ,vq,wg)U]7e-pa^-*V 

9=2, 4,6 

(^33,  iTn,  (723)  =  Z  P{D\q,  D2q,  D3q)V\qe.  pr>dx3  r3) 

<7=2. 4, 6 
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Specialization  of  formal  solutions  to  infinite  media 


At  this  point,  a  formal  solution  of  the  fi<*l<|  equation  in  a  monoelinir  medium 
has  been  presented.  The  amplitudes  / ’,9  are  the  unknowns.  The  amplitudes 
will  now  he  determined  by  implementing  the  artificial  interface  conditions 
(1.14-1.15).  To  this  end.  if  ( 28)  is  subjected  to  t  he  condi t  ions  (11)  and  (29)  to 
(15).  one  finally  solves  the  displacement  amplitudes  as 


c.5  -  t*)F{p)QH'lcnDHmp) 

(30a  i 

C,  -  vs)F(p)Q/('2r^D,imp) 

(30b) 

>’3  -  (’i ) F( p)Q / (2<’:m D ,imp) 

(30c) 

where 


Dum  =  e, ( o3 w-.i  -  oS(fs)  -f  i’.)(o sh;5  -  rv,  re, )  +  t'5(o,  re,  -  n-3  ».’•,)(  2. 31  ) 


It  is  interesting  to  note  that  Dum  =  0  defines  an  equivalent  ChristofFel  char¬ 
acteristic  equation  for  the  propagation  of  bulk  waves  in  the  medium.  With 
these  solutions  for  the  wave  amplitudes,  solutions  in  the  upper  region  can  be 
written  in  terms  of  q  —  1,3.5  as 


(«t.u2,ua)=  £  (-1  ,-vq.wq)i:iqc-^-^ 

9=  1 .3.5 

(<T33.  ^13,  rr-;3  )  =  £  p(-D  ,7.  Dlq .  D3n)l'lqc~ra',{x'~x*) 


(2.32) 


q=  1 .3,5 

In  summary,  solutions  (28)  and  (32)  with  the  amplitude  solutions  (30)  uniquely 
define  the  propagation  fields  in  the  lower  and  upper  “artificial”  semi-spaces. 
In  other  words,  their  combination  constitute  the  total  solutions  for  tin1  infinite 
medium. 


Y  x' 

A3  »  a3 


Figure  2.2:  Semi-infinite  media 

Specialization  of  formal  solutions  to  half-spaces 

One  now  adapts  the  solutions  of  the  infinite  media  (2.28)  and  (2.32)  for  solution 
of  the  case  where  the  free  boundary  is  intercepts  the  propagating  waves  at 
some  arbitrary  location  parallel  to  the  plane  13  =  0-  Assuming  that  the  free 
boundary  is  located  at  x3  =  —d  as  depicted  in  figure  2.  The  implication  is  that 
the  free  boundary  is  located  in  the  upper  region  and  thus  can  only  interfere 
with  the  propagation  fields  in  the  negative  13  direction.  For  this  case,  the 
solution  (2.32)  will  constitute  an  incident  wave  on  free  surface.  As  a  results 
waves  will  reflect  from  the  free  boundary  and  propagate  in  the  positive  x3 
direction.  Thus,  appropriate  formal  solutions  for  the  reflected  waves  can  be 
adapted  from  the  solution  (2.28)  in  accordance  with  (Note  now  that  X3  does 
not  appear  because  solutions  are  referred  to  the  origin  x3  =  0) 

=  £  (1  (2.33) 

<7=1 ,3,5 

(^.^3.^3)  =  £  P(Ou,,D2q,D3q)Ullte-pa‘>Xi  (2.34) 

7=  1,3,5 
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With  this  the  total  solution  for  the  semi-space  (designated  with  superscript 
“s')  which  is  required  to  satisfy  the  stress  free  boundary  condition  is  obtained 
by  superposing  the  incident  waves  and  reflected  waves  in  accordance  with 


(«?,* 3^5)  =  E  (-!•■ 


\q 


q—  1,3,5 


+  £  {1 


9=  1,3,5 


(2.-55) 


<7=  1,3,5 

+  £  p(D,„D!„D3, 
g=  1,3.5 


(2.36 ) 


The  boundary  condition  is  given  by 


*13*  =  <723*  =  £33*  =  0.  at  x3  =  —</ 


(2.37) 


By  imposing  the  boundary  conditions  (2.37)  on  Eq.(6),  linear  system  is  ob¬ 
tained  by 


(D  „  D13  D1S\  /£/[,£,  \  (RA 

D2\  £*23  D25  U[2  E:i  =  /?3 

t  £*31  £*33  £*35  /  \  £[5  /  \  / 

where 

«1  =  +  t/1^13^3  +  u"dVjed 

«3  =  {U$D2,E\  +  6^£>23£3r  +  uSd25ED 
Rs  =  (1/}?D3i£T  +  £','^33  £3  +  U$DMEl) 


(2.38) 


(39a) 

(39b) 

(39c) 


Using  Cramer’s  rule  solve  for  the  reflected  amplitudes  obtaining 

£'"n  —  7;  — (/2iCii  —  R3G2 1  -I-  R^Cri] ) 

L'sm  E\ 

£[3  =  77  —  (/?iG'i3  —  R3G23  +  £5  ££33) 

U sm  ^3 


(40a) 

(40b) 
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(40c) 


ilh  ~  — - —  (/?i6'15  —  /^3^'25  +  R5G3Z,} 

Usm  &S 


Eg  =  ep0lqd,  Er  =  e-p^<r?+!,\  q—  1 .3.5 


(2.41) 


Dsm  =  D\\G\1  —  D'i\  G 2\  +  D31G31 


(2.42) 


Where  the  Gqr  are  given  in  Appendix  B.  Note  that  Dsm  =  0  defines  the 
characteristic  equations  for  the  propagation  of  Rayleigh  (surface)  wave  on  the 
free  surface  [1] 

2.5  Cargniard-DeHooop  contour  variation 

Infinite  media 

Now  consider  the  transformations  back  to  the  time-space  domain  using  Carg- 
niard  DeHoop  method.  Choose  F(t)  as  the  delta  function  S(t)  Laplace  trans¬ 
form  of  which  is  equal  to  1.  Then,  consider  the  Laplace  transform  of  u\.  u\  is 
obtained  by 


=  r  r  Ux  3 

J  —  oo  J~o o 

/oo 

Uise~p^asl3~JVX')dT) 

-OO 


)dj]+ 


(2.43) 


Since,  next  one  wishes  to  carry  out  the  integration  with  respect  to  rj  along  a  cer¬ 
tain  contour  in  the  complex  r)-plane  that  deviates  from  the  positive  imaginary 
axis,  extend  the  definition  of  the  relevant  integral  into  the  complex  //-plane 
by  analytic  continuation  away  from  the  imaginary  7/-axis.  This  deformation 
of  the  integration  contour  is  only  valid  if  no  singularities  are  crossed.  One 
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can  prove  (Burridge  1970)  that  the  only  singularity  in  the  //-plane  are:  (a) 
branch  points  on  the  imaginary  axis  where  up-  and  down-going  waves  meet: 
(b)  possible  branch  points  off  the  imaginary  axis  where  other  two  upgoing 
or  downgoing  waves  meet;  (c)  poles  on  the  imaginary  axis  to  the  lowermost 
point,  and  to  the  uppermost  branch  point  (The  present  of  poles  is  linked  to  the 
existence  of  interface  waves  at  one  or  more  of  the  interfaces  that  are  present 
in  the  configuration;  the  poles  do  not  occur  in  the  case  of  an  infinite  homoge¬ 
neous  solid).  The  singularities  of  type  (b)  are  no  obstruction  to  the  contour 
deformation  process  as  long  as  summations  of  the  contributions  of  the  relevant 
three  generalized  rays  are  taken  into  account  by  adding  the  results  from  the 
three  contour  integrals.  Thus,  in  the  deformation  process  one  only  has  to  as¬ 
sure  that  one  does  not  cross  any  of  the  branch  points  on  the  imaginary  r/-axis. 
Now,  change  the  contour  of  the  //  integration  from  the  positive  imaginary  axis 
in  the  complex  //-plane  to  the  Cargniard-DeHoop  contour,  parameterized  by 
the  real  time-variable  t  running  from  tq ,  the  arrival  time  of  a  certain  wave,  to 
infinity.  This  integral  can  be  recognized  as  the  Laplace  transforms  of  certain 
explicit  functions  of  time,  thus  the  inverse  transforms  by  inspection  as  detailed 
in  Achenbach[32],  The  integration  in  the  complex  //-plane  is  carried  out  along 
the  (in  general,  six  different)  paths  where 

t  =  ax^  —  jrjxi  (2.44) 

with  t.  real  and  positive.  The  paths  that  are  traced  by  these  different  //  as  t 
increases  are  the  Cargniard-DeHoop  contour.  The  Cargniard-DeHoop  method 
is  based  on  the  following  elementary  property  of  the  one-sided  Laplace  trans¬ 
form:  for  given  Laplace  transforms, 

to o 

il,  -  /  ul{xl,x3,t)e~pldt  (2.45) 

Jtq 
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the  inverse  of  the  integral  is  exactly  of  the  form  of 


Ui  =  Ui(xx,x3,t)H{t  -  tq) 


(2.46) 


where  H(t-tq)  is  the  Heavyside  step  function.  Equation  (2.44)  is  written  with 
respect  to  a  by 


t  +  jqx  i 
a  -  - 


(2.47) 


*3 

Substitution  of  (2.46)  into  the  characteristic  equation  (40)  yields  the  sixth 
order  polynomial  defining  the  Fourier  parameter  q 


ft(zi,x3,*,ij)  =  0 


(2.48) 


where 


fl(xi,x3,f,7?)  =  Beqb  +  B3q5  +  B4t)4  +  B3q3  +  B2q2  +  Bxq  +  50(2.49) 


The  coefficient  /?,■  is  given  in  Appendix  A.  The  distinct  six  roots  of  q  are 
obtained  from  (2.47)  with  complex  coefficients,  which  is  the  function  of  the 
spatial  variables  (11,13)  and  time  (t).  The  six  root  of  q  is  composed  of  three 
parabolas  with  respect  to  time  t  for  a  certain  position  (1^x3),  when  they 
are  plotted  as  real  q  versus  imaginary  q,  and  they  are  symmetric  about  the 
imaginary  q  axis.  In  the  interesting  range  of  time,  the  q  curves  form  three 
branches.  Each  of  the  parabolas  is  associated  with  six  distinct  roots  of  a  from 
the  characteristic  equation  (20),  three  of  which  correspond  to  the  lov.er  half¬ 
space  and  others  to  the  upper  half-space,  because  of  the  boundedness  of  the 
waves.  Each  of  the  three  q  represents  a  separate  wavefront. 

Now  the  inverse  Laplace  transform  can  be  obtained  by  mere  inspection  of 
(3.25)-(3.26)  to  be 
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(38a) 
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■('h  HirW  ~  h) 
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where  t j,  t2,  and  t3  are  the  arrival  times  of  the  various  wavefronts.  In  other 
words,  they  are  the  times  that  correspond  to  the  values  of  the  imaginary  ij- 
axis  intercepts  of  the  three  branches  of  r/.  The  notation  r/,+  denotes  the  branch 
of  tj,  to  the  right  side  of  the  imaginary  r/- axis  and  i/~  denotes  the  branch  of 
T)l  to  the  left  side  of  the  imaginary  ?/- axis,  by  using  the  same  technique  the 
displacement  u2  and  u3  are  obtained  by 
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Half-space  media 

Employing  the  above  procedure,  integral  transform  solutions  of  Eq.(35)  in  the 
half-space  are  inversely  transformed  back  to  the  space-time  domain.  However, 
It  is  impossible  to  consider  the  transient  waves  in  the  whole  half-space  by 
using  Cagniard-DeHoop  method,  since  the  order  of  the  exponential  is  the 
linear  function  of  q1,c>3,  and  a3.  Fortunately,  one  can  investigate  the  surface 
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waves  on  the  free  surface,  i.e.,  x:!  =  —d.  The  retlerte<l  waves  of  Eq. for  the 
case  of  x-j  =  0  and  x:J  =  —d  that  is  designated  as  superscript  rs  are  simplified 


—t'33 DumDSJn  * 

— w- — ‘ 


dt'33  DumDsrn  ,  r  s 


2  C33  D  lLr,i  Ds 


Kin  =  (t'5  -  t’3 )(DnC"t  ~  £>2,67  -  D;iiG'l)t~’'“'‘l  + 

(t’t  -  ^)(DUicr;  -  Dncr:  - 

(t’3  -  tti)(D,5G7  -  D^G'2  -  D^GDt-^  (2.50) 
«2m  ~  (f,s  ~  <->;})(  Z)n  6',  +  DnG\  +  D-\\ 63  )c 

(t’l  -  t.’S)(Dl3G7  +  D2367  +  D3;i67)c-^'1',+ 

(t’3  —  t’,)(D|sG,  +  D-i 562  +  D:i 563)6  pr‘4^  (2.ol) 

=  (t’5  -  t’3 ) ( D„67  -  D21C"  -  D3lG%)e-»'"d  + 

(t’l  -  t»5)(D1367  -  D23C2  -  A0G3  )e~p"v<+ 

( t’3  -  t’,)(D1567  -  D,567  -  D35G3)t -p"v<  (2.52) 


c?  =  £  c,r.  g;  =  £  G,,«,,  cr  =  £  c,r», 

9=1  9=1  9  =  1 


(2.53) 


where  G?r  is  given  in  Appendix  C.  Now  the  particle  displacements  are  obtained 
by  superposing  the  incident  waves  on  the  free  surface  of  Eq.(29)  and  reflected 
waves  on  the  free  surface  of  Eqs. (50-52).  Now,  the  remaining  task  is  to  apply 
Cagniard-Dehoop  method  the  same  as  infinite  media.  The  result  for  tz,  are 
given  by 
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The  remaining  displacements  u2  and  us3  will  be  obtained  by  similiar  procedures. 


2.6  Numerical  illustrations  and  discussions 


In  this  section  numerical  illustrations  of  the  above  analysis  are  given.  Let's 
choose  for  our  illustration  a  cubic  material.  The  microstructure  of  this  partic¬ 
ular  material  has  been  extensively  studied  and  models  predicting  its  effective 
anisotropic  properties  are  available  elsewhere  (lj.  The  cubic  material  of  In  As 
is  given  with  respect  to  the  reference  coordinate  system  x'  as 
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which  confirms  the  earlier  conclusion  that  the  transformed  matrix  takes  the 
format  of  monoclinic  symmetry.  Having  chosen  the  material  let’s  summarize 
a  "flow  chart”  like  procedure  for  our  subsequent  calculations. 

After  specifying  the  azimuthal  angle  $  (namely,  the  line  load  direction) 
we  proceed  in  the  first  step  to  evaluate  r/’s  from  the  sixth  order  polynomial 
(49)  for  given  time  and  location  (/|,i3),  While  time  (t)  is  n. arching  front  0. 
rf s  have  three  pairs  of  complex  conjugates  for  t  >  two  pairs  of  complex 
conjugates  for  t2  <  t  <  1 3,  and  a  pair  of  complex  conjugates  for  <  t  <  t,. 
tq  represents  the  arrival  times  of  three  kinds  of  wave.  As  the  tj  doesn’t  have 
a  pair  of  conjugate,  displacements  and  stresses  are  equal  to  zero.  The  arrival 
times  of  three  kinds  of  waves  are  collected  with  respect  to  arbitrary  direction 
of  wave  propagation  in  xt  —  x3  plane.  The  arrival  times,  the  energy  slowness, 
corresponding  to  arbitrary  angle  of  <p  for  the  cubic  material  are  shown  in  figure 
5-8.  However,  the  cubic  material  has  cuspidals  in  the  plane  of  arrival  time 
that  represents  lacunae,  in  which  the  elastic  energy  is  zero.  In  the  cuspidal, 
the  material  has  as  many  as  five  arrival  times.  The  number  of  arrival  times, 
therefore,  are  dependent  on  materials.  This  agrees  with  [24].  The  inverse  of 
the  plane  of  arrival  time  is  the  plane  of  the  elastic  energy,  which  shows  the 
energy  fronts  of  the  wave  propagation  in  the  infinite  media.  Figure  5-8  shows 
the  variation  of  energy  front  corresponding  to  azimuthal  angle.  This  energy 
plane  is  exactly  same  as  that  calculated  by  the  plane  of  slowness. 

Next  step  is  to  solve  the  characteristic  equation.  One  then  has  to  decide 
which  of  the  a,,  correspond  to  r/, ,  i  —  1,3,5.  If  one  determines  the  a„  from 
ads  that  match  with  a  of  equation  (27),  the  order  of  other  two  o,7.i  7^  j 
does  not  affect  to  the  numerical  results.  The  sign  of  Re(ou)  follows  sign  of 
Re(or„),  However,  if  the  identification  of  the  o’s  is  in  error,  the  final  results 
of  the  displacements  will  clearly  indicate  that  something  is  amiss.  Lacking 


the  proper  identification  of  the  a‘s  the  calculated  displacements  will  contain 
imaginary  parts,  a  violation  the  condition  for  the  existence  of  the  Laplace 
transform.  Therefore,  only  the  proper  identification  of  the  a’s  will  yield  purely 
real  displacements  so  that  the  existence  condition  of  the  Laplace  transform  is 
satisfied.  Of  course  once  the  existence  condition  is  satisfied  then  the  validity 
of  the  identification  scheme  was  required  since  all  of  the  explicit  solutions  for 
cubic  equations  '’mixed”  roots  for  the  various  of  the  parameter  .  Then  the 
displacement  ratios  are  calculated.  Finally,  displacements  are  determined  by 
explicit  solution  of  (31). 

Response  of  infinite  media 

Numerical  results  for  cubic  medum  of  InAs  with  (j>  =  30°  are  shown  in  figure  3- 
4.  Figure  3a  shows  the  slowness  curves  and  the  direction  of  wave  propagation 
10°  and  5°.  Imaginary  part  of  rj  variation  with  respect  to  time  is  shown  in 
figure  3b  or  4a.  Thus  designating  ti,t2,t3,t4,  and  ts  as  the  arrival  times  of 
various  wave  forms  in  order  from  figure  figure  3b  or  4a.  The  wave  along  the 
paths  of  T)  corresponding  to  <  t  <  t3  or  t  >  t3  represent  quasi-longitudinal 
wave.  The  wave  along  a  rj  of  t  >  t4  represent  quasi-vertical  shear  wave.  The 
wave  along  a  r]  of  t  >  f5  represent  quasi-horizontal  shear  wave.  The  time 
interval  of  t2  <  t  <  t3  represents  ‘lacunae’  that  all  the  disturbances  are  zero. 

Figures  5-8  show  numerical  illstrations  corresponding  to  the  azimuthal 
angle  of  the  chosen  several  material  systems  :  (a)  variations  of  the  energy 
flows,  (b)  energy  slowness  curves,  and  (c)  displacement  fields.  In  fact,  the 
energy  slowness  curves  show  the  arrival  times  for  the  various  wave  forms. 
Both  energy  (a)  and  energy  slowness  (b)  curves  are  based  on  the  slowness 
curves  [7],  (c)  present  snap  shots  of  absolute  value  of  radical  displacement 
field  at  fixed  time  (t  =  0.2  microsec).  A  spatial  grid  of  100x100  points  are 


generated  for  the  first  quadrant.  The  remaining  quadrants  is  then  generated 
by  the  mirror  of  the  first  quadrant  .  The  vertical  line  load  is  located  at  the 
origin.  The  origin  is  located  at  the  center  of  the  picture.  Since  the  medium 
is  homogeneous,  the  wave  field  does  not  change  as  it  propagates.  The  color 
scheme  runs  from  white  (minimum}  to  black  (maximum).  In  this  picture  we 
can  clearly  recognize  the  wave  curves  (three  wave  fronts  and  lacunae).  This 
displacement  intensity  pictures  agree  with  the  energy  flows  (a). 

The  variations  of  the  displacement  fields  corresponding  to  azimuthal  angle 
for  the  cubic  material  of  InAs  are  shown  through  figure  5-8.  Very  sharp  dis¬ 
placement  intensity  and  lacunae  is  presented  in  figure  5  in  the  cse  of  azimuthal 
angle  <f>  —  0°.  As  the  azimuthal  angle  is  increased  by  10°,  the  third  wave  front, 
quasihorizontal  shear  wave  is  more  clear  and  the  lacunae  disappears.  We  can 
not  get  the  good  numerical  results  arround  90°  direction  because  of  contour 
sigularity  of  integral  (47). 


Response  of  half-space  media 


In  this  section,  numerical  results  in  half-space  media  are  discussed.  The  first 
case  is  to  evaluate  the  response  generated  by  an  internal  line  load.  One  choos¬ 
ing  for  our  illustration  a  cubic  material  of  the  perturbed  steel,  the  material 
properties  are  given  with  respect  to  the  reference  coordinate  system  x't  as 
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with  p  —  1.7  g/cm3.  Figure  9  show  a  slowness  curves  (a)  and  displacements 
versus  time,  iq  (b)  and  u3  (c),  for  4>  =  30°.  The  receiver  is  located  at  (2.0, 
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Figure  2.6:  Displacement  field  due  to  a  line  load  in  the  infinite  cubic  medium 
of  InAs  (<f>  =  10°)  :  (a)  Slowness  surface  ;  (b)  Energy  slowness  ;  (c)  Energy 
flow  ;  (d)  Snap  shot  of  displacement  field  at  t=0.2  microsec 


Figure  2.7:  Displacement  field  due  to  a  line  load  in  the  infinite  cubic  medium 
of  In  As  ((f)  =  20°)  :  (a)  Slowness  surface  ;  (b)  Energy  slowness  ;  (c)  Energy 
flow  ;  (d)  Snap  shot  of  displacement  field  at  t=0.2  microsec 


0.0,  0.2)(mm).  Sharp  variations  of  displacements  are  due  to  the  arrival  of 
the  various  waves.  As  time  increases,  the  displacements  approach  to  zero, 
which  is  consequence  of  the  direc  delta  function.  The  res>‘!,s  in  figure  9  are  as 
expected  physically.  A  strong  peak  is  observed  that  shows  the  existence  of  a 
surface  or  a  Rayleigh  wave.  We  can  clearly  recognize  quasi-longitudinal  wave. 
2  quasi-shear  waves,  and  a  surface  wave. 

The  second  case  is  to  numerically  calculate  the  response  generated  by  the 
line  load  on  the  surface.  The  surface  is  free  from  tractions  except  the  line 
load  along  the  x2-axis.  Many  researchers  (Taylor,  Krout,  etc)  investigated 
this  problem  for  the  various  line  load  such  as  direc  delta  function  or  heavyside 
step  function  with  respect  to  time.  Analytical  results  for  the  heavyside  step 
function  surface  line  load  in  the  half-space  media  are  given  in  [29].  We  present 
numerical  illustrations  of  the  snap  shot  for  the  displacement  fields  generated  by 
direc  delta  line  load.  Figure  10  present  snap  shots  of  radical  displacement  field 
at  fixed  time  ( t  =  0.2)  for  isotropic  steel  (a),  perturbed  steel  of  <p  =  0°  (b),  and 
perturbed  steel  of  0  =  30°  (c).  A  spatial  grid  of  100x100  points  are  generated 
for  the  first  quadrant.  The  fourth  quadrant  is  then  generated  by  the  mirror 
of  the  first  quadrant.  The  color  scheme  runs  from  white  (maximum)  to  black 
(minimum).  In  these  pictures  we  can  clearly  recognize  the  wave  curves  (head 
wave,  surface  wave,  and  three  bulk  wave  forms).  The  surface  wave  attribution 
is  shown  at  inside  of  longitudinal  wave  front  and  near  to  the  horizontal  plane 
surface. 
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Figure  2.10:  Snap  shot  of  displacement  field  due  to  a  surface  line  load  of  the 
semi-  infinite  medium  at  t=0.2  microsec  :  (a)  Displacement  field  for  steel  ;  (b) 
Displacement  field  fur  perturbed  steel  (  f=1.5,  d  =  0°)  ;  (c)  Displacement  field 
for  perturbed  steel  (f=1.5,  <i>  —  30°) 
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Appendix  A 

Various  coefficients  of  characteristic  eq.(18)  are  given  by 

An  =  c55a2  -  cn7)2  -  p 
A12  -  c45a2  -  Ci6??2 
Aia  =  ~jr]oc{c  13  +  C55) 

A22  =  c44a2  -  c^T}2  -  p 
A23  =  —jvot(c3&  +  c.  45) 

A .33  =  c33a2  -  C53172  -  p  (2.55) 


71 


Ax  =  Dx 

Ay  =  Dll)2  +  Dip 

j43  =  D^rj4  +  Dapr)~  +  D$p~ 

A4  =  D?*?6  +  Dspp4  +  Dgp‘7)2  +  p' 3  (2.56) 

Appendix  B 

Various  coefficients  of  polynomial  for  Cargniard-DeHoop  contour  are  given  by 

B6  =  D2x\x\  —  D4X3X1  —  +  £)7x3 

B5  =  Cxjt 

B4  =  C2t 2  +  Dix\x\  —  D5x*x\  +  D8x  3 
B3  —  C3]t3  +  C^jt 
B2  -  Cot  '  +  C6t2  -  £>6X3X1  +  £>9x3 
Bx  =  6Dxjxxt5  +  4Z)3jx3Xit3  +  2D6x8xxt 

B0  =  Dxt6  +  DixlC  +  £>6^3  +  P3x63  (2.57) 

Cx  =  6DiXi  +  2£>4x3X]  —  4£>2a:3Xj 
C2  =  lSD^  —  6£>2x3Xi  +  £>4x3 
C3  =  4D2x3Xi  —  20D4XJ 
C4  =  — 4D3x3Xi  +  2D5X3X1 
C5  =  D2xl-\$Dxx\ 

C6  =  —6  03x3x^  +  D.5X3  (2.58) 

DX  ~  C33C45  —  C33C44C55 

D2  —  C\  ]  033044  +  C55 F\  —  Ci6C33C45  —  C45F2  +  (ci3  +  C,55)F3 

Dn  =  [055(033  +  C44)  +  C33C44  —  c24S]p 

Da  =  [(Ci3  +  Caa)F.\  +  C45C16O55  +  C16F2  —  C55C66  —  CnFi] 

D$  —  [(c13  +  C55)2  +  2C45C1S  —  Fx  —  055(055  -f  cee)  —  Cj  1  ( C33  +  o4.})]p 
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